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Abstract 

The effects of nonuniversality of gaugino masses on dark matter are ex- 
amined within supersymmetric grand unification, and in string and D brane 
models with R parity invariance. In SU(5) unified models nonuniversality 
in the gaugino sector can be generated via the gauge kinetic energy function 
which may depend on the 24, 75 and 200 dimensional Higgs representations. 
We also consider string models which allow for nonuniversality of gaugino 
masses and D brane models where nonuniversality arises from embeddings 
of the Standard Model gauge group on five branes and nine branes. It is 
found that with gaugino mass nonuniversality the range of the LSP mass can 
be extended much beyond the range allowed in the universal SUGRA case, 
up to about 600 GeV even without coannihilation effects in some regions of 
the parameter space. The effects of coannihilation are not considered and 
inclusion of these effects may further increase the allowed neutralino mass 
range. Similarly with the inclusion of gaugino mass nonuniversality, the 
neutralino-proton (x — p) cross-section can increase by as much as a factor 
of 10 in some of regions of the parameter space. An analysis of the uncer- 
tainties in the quark density content of the nucleon is given and their effects 
on x ~ P cross-section are discussed. The predictions of our analysis includ- 
ing nonuniversality is compared with the current limits from dark matter 
detectors and implications for future dark matter searches are discussed. 



1 Introduction 

In this paper we study the effects of nonuniversality of gaugino masses on dark 
matter in SUGRA models Q and in string and D brane models. Nonuniversal 
gaugino masses arise quite naturally in supergravity and string unified theories. 
Thus, in N=l supergravity the kinetic energy and the mass terms for the gauge 
fields and the gauginos are given by|l| 
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Here X a are the gaugino fields, G = — lnf/t 6 !^!^*] — K 2 d, where W is the superpo- 
tential, d(z, z*) is the Kahler potential where z a are the complex scalar fields, and 
k = (8ttGn)~^ = 0.41 x 1CT 18 GeV -1 where Gn is Newton's constant. The gauge 
kinetic energy function f a p in general has a non-trivial field dependence involving 
fields which transform as a singlet or a non-singlet irreducible representation of the 
underlying gauge group. After the spontaneous breaking of the gauge symmetry at 
the unification scale to the Standard Model gauge group SU (2) L x U{1) x SU (3)c, 
one needs to carry out a rescaling of the gauge kinetic energy in the sector of 
the gauge group that is preserved. This rescaling generates a splitting of the 
SU{2) L x U{1) x SU(3)c gauge couplings at the unification scale M x and one 
has!, |, 1 

a~\M x ) = a x \M x ) + £co r < (2) 

r 

where n[ characterize the Higgs vacuum structure of the irreducible representation 
r and Cq t parametrize its relative strength. After rescaling the gaugino mass matrix 
takes the form 

map = \e G/2 G a (G- l ) b a (df: 7 /dz* b )f^ (3) 

Here one finds that the contribution to nonuniversality of the gaugino masses is 
controlled not only by the nature of the GUT sector but also by the nature of 
the hidden sector. Because of this the splitting of the gaugino masses at M x is 
characterized by a set of parameters c r different from cq t . Thus after the breaking 
of the unified gauge group we parametrize the gaugino masses at M x by m^O) 
where §, |] 

rhi(0) =mi(l + ^c r <) (4) 

r 

The effect of Cq t on the gauge coupling unification has been discussed in the pre- 
vious literature J§] and we do not discuss it here. In the analysis of this paper we 
assume unification of gauge couplings at M x and the refinement of including co r 
correction will not have any significant effect on our analysis. 

In SU(5) the gauge kinetic energy function f a p transforms as the symmetric 
product of 24 x 24 and contains the following representations 

(24 x 24) symm = 1 + 24 + 75 + 200 (5) 

The singlet leads to universality of the gaugino masses while the non-singlet terms 
will generate nonuniversality. We consider models where we have a linear combina- 
tion of the singlet and a non-singlet representation, i.e., 1 + 24, 1 + 75 or 1 + 200. 
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The quantities n\ for the representations 1, 24, 75, 200 are listed in Tablel|| [|. [|. 
Some phenomenological aspects of nonuniversality of gaugino masses have recently 
been discussed p, [7|, ||. We focus here on their effects on event rates in the direct 
detection of dark matter (see Ref.[[| for previous work on the effects of gaugino 
mass nonuniversality on dark matter). 

Techniques for computing the event rate in the scattering of neutralinos off 
nuclear targets has been discussed by many authors[ID|. We follow here the proce- 
dures discussed in Ref.|ITJ. In our analysis we impose the h — > s + 7 constraint ]T2| 
and the bounds on SUSY particles from the Tevatron and LEP[13|. Specifically 
we take for the lower limits m. 



> 94 GeV, m x > 33 GeV, m fR > 71 GeV, 



171; 



> 87 GeV, m g > 190 GeV, m h > 113.5 GeVfTj and for the b -> s + 7 



i.e., we 



branching ratio we take a 2a range around the current experiment fll5 
take 2 x 10~ 4 < B(b — > s + 7) < 4.5 x 1CT 4 . The quantity that constrains 
theory is Q x h 2 where Q x = p x j p c , where p x is the neutralino relic density and 
p c = 3Hq/8t[Gn is the critical matter density, and h is the value of the Hubble 
parameter Ho in units of 100 km/sMpc. The current limit on h from the Hubble 
Space Telescope is h — 0.71 ± 0.03 ± 0.07|16| and recent analyses of fl m give|l7] 



Vl m = 0.3 ± 0.08. Assuming Vl B ~ 0.05, one gets 



Q v h 2 = 0.126 ±0.043 



(6) 



In this analysis we make a a somewhat liberal choice for the error corridor on Q x h 2 , 
i.e., we choose 0.02 < Q x h 2 < 0.3. The choice of a more restricted corridor does 
not affect the general conclusions arrived at in this analysis. In the theoretical 
computation of the relic density we use 

rl/2 
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where j is the reheating factor, Nf is the number of degrees of freedom at the 
freeze-out temperature Tf and Xf = kTf/m x . In determining J{xf) we use the 
method developed in Ref.jTjJ. The role of J in the context of nonuniversalities 
will be elucidated in Sec. 5. A number of effects on neutralino dark matter have 
already been studied. These include the effects of nonuniversality of the scalar 
masses at the unification scale [19|, |20|, effects of variations of the WIMP velocity [pT 



22 , p3fl , effects of rotation of the galaxy [[24]], effects of CP violation with EDM 
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constraints [|25|], and effects of coannihilation [ 26 ] . The focus of this analysis is 



on the effects of nonuniversality of the gaugino masses on dark matter. In the 
present analysis we do not include the effects of coannihilation. These effects 
become important when the NLSP mass mi lies close to the LSP mass, i.e., Aj = 
{rrii I 'm x —l) < 0.1. We have identified several regions of the parameter space where 
coannihilations involving fx, §r, the light chargino xt '■> or the next to the lightest 
neutralino x° occur. However, as we stated above we do not consider coannihilation 
in this paper and thus eliminate such regions of the parameter space by imposing 
the constraint Aj > 0.1. An analysis in this region requires a separate treatment 



and will be reported elsewhere p7| . Recent analyses[gg, g{| have pointed to the 
uncertainties in the quark masses and in the quark content of the nucleon that 
enter in analyses of dark matter. We give in this paper an independent analysis 
of the errors in the quark densities and compute their effects on dark matter. 

The outline of the rest of the paper is as follows: In Sec. 2 we discuss the basic 
formulae used to compute the neutralino-proton cross-section. An analysis of errors 
in the quark densities that enter in the scalar a xv cross-section is also given. In 
Sec. 3 we first give an analysis of a xp (scalar) for the universal SUGRA case and 
analyze the effect of errors on the quark densities of the proton on it. We then 
discuss three nonuniversal scenarios where we consider admixtures of the singlet 
with the 24 plet, the 75 plet and the 200 plet representations for the gauge kinetic 
energy function. In Sec. 4 we extend the analysis of a xp (scalar) to the case of the 
O-II string model and a brane model based on 9 branes and 5j branes. In Sec. 5 we 
discuss the origin of the enlargement of the allowed LSP domain consistent with 
the relic density constraints due to the presence of nonuniversalities. Conclusions 
are given in Sec. 6. In Appendix A we give an analytic solution of the sparticle 
masses using the one loop renormalization group equations including the effect 
of gaugino mass nonuniversalities. Using results of Appendix A we compute in 
Appendix B the effects of nonuniversalities on the fi parameter. The analytic 
results of Appendices A and B provide a deeper understanding of the gaugino- 
mass nonuniversality effects discussed in Sees. 3, 4 and 5. 



2 Neutralino-proton cross-section 

For heavy target nuclei such as germanium the neutralino-nucleus scattering cross- 
section is dominated by the scalar part of the neutralino-quark interaction and it 
is the quantity a xp (scalar) on which constraints have been exhibited in the recent 
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experimental works ^T|. For this reason we focus in this paper on the analysis 



of a xp (scalar) . The basic interaction governing the x~P scattering is the effective 
four-fermi interaction given byp2 



A// = XI^XIY^Pl + BP R )q + Cxxm q qq + Dxj 5 Xm q q-f 5 q + Ex^X^qq 

+Fxxm g qi-f5q (8) 

where the interaction relevant to our analysis is parametrized by C . The x ~ P 
cross-section arising from scalar interactions is given by 

Au 2 2 
a xp (scalar) = £ f?Q + -(1 - £ /f) £ C a f (9) 

i=u,d,s i=u,d,s a=c,b,t 

Here // r is the reduced mass, ff (i=u,d,s quarks) are defined by 

m pfi =< P\m qi qiqi\p > (10) 

and C is given by 

C = C h0 + C H o + C f ~ (11) 

where C^o, C#o are the contributions from the s-channel h° and H° exchanges and 
Cj is the contribution from the t-channel sfermion exchange. They are given by 



AMy/M^o sm p [cos p) 



~ f ,\ g 2 sin a(cosa) . . 

4My^ilf^ sm p(cosp) 



C / ( "' d) = "i^ M^f " i^M^Mf (14) 

Here (u,d) refer to the quark flavor, a is the Higgs mixing angle, and C q L, C' qL etc. 
are as defined in Ref. E2J, and a and p are defined by 



o = X^X^q — tan 9wX* ) cos a + X^X^q — tan 0^X* O ) sin a (15) 
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p = — X| (Xj — tan 9wX{q) sin a + X^X^ — tan 9wX* Q ) cos a (16) 
where X n o are the components of the LSP 

X = X* B + X; W 3 + X^Ht + Xl H 2 (17) 

We discuss now the amount of uncertainty connected with the determination of 
ff . The quantities that are used as inputs are cf-kN, x , and £ defined by 

< p\2~ l {m u + m d ){uu + dd\p >= cr n N, (18) 

a < p\uu + dd — 2ss\p > , . 

x = = n= T7i ( 19 ) 

(?ttN < p\uu + dd\p > 

and 

< p\uu + > 
We can determine ff in terms of these and find 

tp 

■1 a 



r: 





(1 + 




m u + m d 


m p 


m d 


(i-0 




m u + m d 


m p 


m. 


l-x) 






m u + m d 


m p 



fl = —(l-x)— (21) 

m u + m d m p 

A similar analysis holds for the neutralino-neutron scattering and one can de- 
termine /" in terms of £, x, a^N as follows 

f" 

J u 

en 
Id 





(i-0 


0~ttN 


m u + m d 


m p 


m d 


(i + 


&ttN 


m u + m d 


m p 


m s ^ 


l-x) 




m u + m d 


m p 



ff= 8 (l-x)— (22) 
We note in passing that from Eqs.21 and 22 one has the relation 

flfl = fU2 (23) 

which holds independent of the details of the input parameters. We discuss now 
the numerical evaluation of ff and ff. The various determinations of ovjv, o"o and 
x using analyses of Ref.[B3|, [34], |35L [361 [37], |28| are summarized in Table2. For a n N 
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one finds an average value of 48 ± 9 MeV, and for do an average value of 35.5 ± 6 
MeV. These give x = 0.74±0.25. Further, there are two independent lattice gauge 



determinations^, |37j of y = l — x which we list in Table2. In recording the result 



of x for Ref. 137] we have reduced the y value by 35% as discussed in Ref. p6 



The average of these lattice gauge calculations gives x = 0.61 ± 0.08. Taking the 
average yet again of this x and of ao/cr-KN we get the average x listed in Table2. 
In addition to the above we need to determine the symmetry breaking parameter 
£. Here as in the work of Ref.|29| we use the analysis of Ref.|38[] on baryon mass 
splittings to obtain 

, = (5- + 5°-s + -s-)* (24) 

s S- + 2° + £+ + E- - 2m p - 2m n v ; 

where x is as defined by Eq.(19). Numerically one finds £ = 0.196a; and on using 
Table 2 we find 

f = 0.132 ±0.035 (25) 

In addition to the above one needs the ratios of the quark masses for which we 
use 



34 



Tfl 777, 

— = 0.553 ±0.043, — = 18.9 ±0.8 (26) 
m d m d 



On using Eqs.(21),(25),(26), and Table 2 we find 

fl = 0.021 ± 0.004 

CP 



Similarly for /" we find 



f d = 0.029 ± 0.006 
ff = 0.21 ±0.12 (27) 



= 0.016 ± 0.003 
fl = 0.037 ± 0.007 

/ s n = 0.21 ±0.12 (28) 

For the more general case of neutralino-Nucleus (x — N) scattering one has 

Am 2 2 
a x _ N (scalar) = —L(Z £ f?C t + -Z(l- £ /f) £ C a 

i=u,d,s i=u,d,s a=c,b,t 

+(A-Z) £ ff Ci + l(A-Z)(l- Yl fn E C a ) 2 (29) 

i=u,d,s i=u,d,s a=c,b,t 

Using Eqs.21 and 22 we write the above in the form 
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Am 2 r A £ m u C u + m d C d £ .m u C u - m d C d 

a x - N (scalar) = (/ ■ + f£A ■ + fC s 

ix m u + m d m u + m d 

where A = (2Z — A) /A, f — f s , f — cr^/m p and numerically on using Table 2 
we have 

/ = 0.05 ±0.01 (31) 

We note that while the \— p and x~ n cross-sections depend on £, the £ dependent 
term has a cancellation in (x — N) cross-section because of the (2Z — A) factor. 
Further, if the target nucleus has A = 2Z, i.e., A = 0, then the £ dependent term 
will drop out of the x ~ N cross-section. Because of the above it is experimentally 
better to plot the % — N cross-section rather than the x ~ P cross-section as is 



currently the practice 30, 31 . 



3 Dark matter in GUT models with gaugino mass 
nonuniversality 

As discussed at the beginning of this section we consider here models where the 
nonuniversalities arise from admixtures of the singlet with the 24 plet, the 75 plet 
and the 200 plet representations. However, we begin first by exhibiting the result 
for the universal SUGRA case. The soft SUSY breaking sector of the theory, under 
the assumption that SUSY breaking is communicated from the hidden to the visible 
sector by gravitational interactions, is parametrized in this case by the universal 
scalar mass m , the universal gaugino mass mi, the universal trilinear coupling A Q 

2 

all taken at the GUT scale, and tan/5 =< H 2 > / < Hi > where H 2 gives mass to 
the up quark and H\ gives mass to the down quark. Throughout this analysis we 
assume that the Higgs mixing parameter \x (which appears in the superpotential 
as /j,HiH 2 ) is determined via the electro- weak symmetry breaking constraint. The 
range of the parameters are limited by a naturalness constraint. We mean this to 
imply that m Q ,mg < 1 TeV, where m~ g is the gluino mass, tan/3 < 25, and A , 
or equivalently A t , the value of A at the electro- weak scale in the top channel, is 
limited by the electro- weak symmetry breaking constraint. For the analysis here 
we choose \i > while for the other \i sign the allowed parameter space for dark 



S 



matter is strongly limited due to the 5 ^ s+7 constraint [12]| . The results for < 
look qualitatively different in that the cross-sections are significantly smaller. 

In Fig.l we plot the maximum and the minimum of a xp (scalar) as a function of 
m x where the parameters are allowed to vary over the naturalness range discussed 
above. The analysis is done for three sets of ff values corresponding to the corridor 
given by Eq.(27). They correspond to = 0.025, = 0.035, ff = 0.33, 

(II) fP = 0.021, /£ = 0.029, ff = 0.21, (III)# = 0.017, f d = 0.023, ff = 0.09. From 
Fig.l we see that the different sets can lead to a variation in a xp (scalar) up to a 
factor of about 5. For the rest of the analysis in this paper we use set (II). 

Next we consider the 1 + 24 model. We find that in this case a xp (scalar) 
typically decreases for positive values of c 24 and increases with negative values of 
C24. This behavior arises primarily from the dependence of the gaugino-Higgsino 
components X nQ of the LSP on the gaugino mass nonuniversality. Thus X n0 are 
sensitive to the gaugino nonuniversality through their dependence on fhi, m 2 and 
\i. In Fig. 2 we exhibit the dependence of X n0 on C 24 for some typical input 
values. The quantity a xp (scalar) depends on the direct product of the gaugino and 
Higgsino components of \- Specifically, a xp (scalar) vanishes if x is a pure Bino. 
From Fig. 2 we see that negative values of C24 increase the Higgsino components 
and hence increase the neutralino-quark scattering and lead to an enhancement of 
a xp (scalar) while the opposite situation is realized for positive values of C24. This 
is what is found in Fig. 3 where we give a plot of the maximum and the minimum 
of a xp (scalar) for the cases c 24 = —0.08, c 24 = and c 2 4 = 0.1 when the soft SUSY 
breaking parameters are varied over the assumed naturalness range as in Fig.l. 
The analysis shows that for c 2 4 = —0.1 an enhancement of a xp (scalar) by as much 
as a factor of 5 can occur as a result of the gaugino mass nonuniversality and the 
allowed range of m x is also increased in this case beyond the range allowed in the 
universal SUGRA case. 

In Fig.4 we give an analysis of the maximum and the minimum of a xp (scalar) 
for the 1 + 75 case for three different values of C75, i.e., C75 = —0.06, c 75 = 
and C75 = 0.04. As for the 1 + 24 case, a xp (scalar) typically increases for negative 
values of C75 and decreases for positive values of C75. Again this can be understood 
by analysing the gaugino-Higgsino components of the LSP as a function of C75. 
Thus here as in the 1 + 24 case one finds that the Higgsino components of the LSP 
increase as C75 decreases and decrease as C75 increases. Because of this there is an 
enhancement of a xp (scalar) for c 75 < 0. A comparison of c 75 = and C75 = —0.06 
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cases in Fig. 4 shows that an enhancement of a XP (scalar) up to a factor of 5 or 
more occurs in this case. As in the case of 1 + 24 here also one finds that the 
allowed range of m x consistent with the constraints is extended beyond the values 
allowed in the universal SUGRA case. 

In Fig. 5 we give an analysis of the maximum and the minimum of the a xp (scalar) 
for the 1 + 200 case when c 20 o = —0.08, c 20 o = and c 20 o = 0.1. In this case the 
dependence of a xp (scalar) on c 2 oo is opposite to that one has in the previous cases, 
i.e., the 1 + 24 case and the 1 + 75 case. Here it is for positive values of c 2 oo 
that a xp (scalar) increases and it is for negative values of c 2 oo that a xp (scalar) de- 
creases. The origin of this reversal lies in n\ (see Table 1) and can be understood 
from the dependence of X n0 on c 20 o- The above implies that the Bino component 
of x decreases and the Higgsino components increase for c 2 oo > while the op- 
posite situation occurs when c 2 oo < 0. This dependence implies that a xp (scalar) 
should increase for c 2 oo > and decrease for c 2 oo < which is what is observed 
in Fig. 5. In this case one finds that a xp (scalar) can increase up to a factor of 10 
or more because of nonuniversality. The analysis also shows that as in the case of 
1 + 24 and 1 + 75 the allowed range of the LSP mass is extended beyond what is 
allowed in the universal SUGRA case. 

One can gain a deeper understanding of the dependence of X n0 on c r and hence 
a deeper understanding of the dependence of a xp (scalar) on c r from studying the 
dependence of the Higgs mixing parameter /i on c r . To appreciate why /i is such 
an important parameter in this discussion it is useful to look at large /x, i.e., the 
case fi 2 » Mf . In this limit one finds that the LSP eigenvector is given by 

X 10 ~ 1 - (Mf/2// 2 )sin 2 w , A 20 ~ -(Mf/2(m 2 - mi)//) sin 26 w sin 2(3, 

X 30 ~ (M z /fi)sin6 w sin(3, X 40 ~ — (M z /n)sin6 w cos/3, (32) 

Eq.32 shows that in the large \i limit x is mostly a Bino and the corrections to the 
pure Bino limit are proportional to (Mf//i 2 ) while the Higgsino components are 
proportional to (Mz/fJ,)- As already pointed out the a xp (scalar) depends on the 
interference of the gaugino and Higgsino components of the LSP, i.e., X i0 x X a0 
(i=l,2; a = 3,4). Clearly then as \fj,\ increases we go more deeply into the pure 
Bino region reducing a xp (scalar). Likewise as decreases x develops larger 
Higgsino components X a0 (a = 3,4), even though it is still dominantly a Bino, 
a xp (scalar) decreases. Thus to gain an insight on the effect of c r on X n0 and 
hence on the effect of c r on a XP (scalar) we need to understand how c r affects fi. 
We address this topic below. 
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In SUGRA models /x 2 is determined via the breaking of the electro-weak sym- 
metry and thus depends on the gaugino mass nonuniversality through the Higgs 
mass parameters (see Eqs.39 and 43 in Appendix A and Eq.50 in Appendix B). 
We can understand the effect of nonuniversality on /x analytically by expanding 
/x for the nonuniversal case around the universal value using c r as an expansion 
parameter 

? =^ + E^ + 0{cl) (33) 

where /io is the value of /x for the universal case. Using the analysis of Appendix 
B the pattern of breaking in the Higgs structure of 24, 75, 200 shows that 

^li >0 ,^4>0,^<0 (34) 
dc 24 oc 75 dc 20 o 

Thus in the neighborhood of c r = a negative value of C24 gives a smaller value 
of |/x| leading to larger Higgsino components in \ and hence a larger a xp (scalar) 
as is observed in the numerical analysis. A similar situation holds for the nonuni- 
versality effects from C75. However, for the nonuniversality effects from C200 an 
opposite situation holds because of the opposite sign of the derivative term as 
given by Eq.34. More generally one finds the same behavior in a larger c r domain, 
i.e., /x 2 4 < /x 2 , for c 2 4 < 0, /x 2 5 < /Xq for c 75 < 0, and nl 00 < /x 2 , for c 20 o > and 
observations similar to those valid for small c r also apply here. These results imply 
that gaugino nonuniversality which makes /x small produces a deviation of the LSP 
from the approximate Bino limit in a direction which leads to a larger value of 
a xp (scalar) . 

Nonuniversality of the gaugino masses also has implications for naturalness. 
One convenient definition of naturalness is via the fine tuning parameter $ intro- 
duced in Ref. [|39| defined by $ = | + j^-. Using this definition one can easily 
compare values of fine tuning for the universal and nonuniversal cases. One finds 
$24 < $0 (c 24 < 0), $ 75 < $0 (c 75 < 0), $ 20 o < $0 (c 20 o > 0), where $24 is 
$ for the case 1 + 24 etc, and $o is & f° r the universal case. Since the correction 
to xx 2 is negative for the case of larger Higgsino components one finds that the 
deviation from the approximate Bino limit is in the direction of a smaller value of 
$ and towards the direction of greater naturalness relative to the universal case. 
Thus a smaller /x leads to a larger a xp and a larger detection rate and also makes 
the model more natural by making $ small. In this sense the more natural the 
SUSY model the larger is the detection rate. 
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4 Dark matter in string/brane models 



One of the main hurdles in the analysis of SUSY phenomenology based on string 
models is that there is as yet not a full understanding of the breaking of supersym- 
metry here. However, there do exist efficient ways to parametrize SUSY breaking 
and one such parametrization is HO 



F s = y/3m^(S + S*) sin^e"* 75 , F i — x/3m| (T + T*) cos 00je -i7 * (35) 

where F s is the dilation VEV, F % are the moduli VEVs, 9 (Oj) parametrizes the 
Goldstino direction in the S (Tj) field space and 75 (7$) is the phase. The 0, obey 
the constraint 01 + 02 + 03 = 1. We begin with an example of the O-II string 
model with the soft SUSY breaking sector parametrized by|}40f 



m 8 = v^m^sinfle - ** - 7 i ecos0e -<ar ) 



m r 



e'(-5 GS )ml, A = -V3rri3 sin 6e~ ias (36) 
2 2 



Here 71 = — P + Sgs, I2 = —1 + 8gs, 73 = 3 + 5gs where 5gs is the Greene- Schwarz 
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parameter which is fixed by the constraint of anomaly cancellation in a given 
orbifold model. Further, as in the GUT analyses we treat fi as an independent 
parameter. The phenomenology of this model has been discussed in Ref. and 
the EDM constraints in Refs. [f42"l, 43]. However, in the analysis below we do not 
impose the accelerator constraints jy] and as in the analysis for GUT models we 
do not assume CP violation and thus set the CP phases to zero. 

The stucture of the soft SUSY parameters for the model discussed above shows 
that the nonuniversality of the gaugino masses is controlled by several parameters 
in this case: 5gs, & and e which play a role similar to the role played by the param- 
eters c r in the case of the GUT models. The presence of several parameters leads 
to many different possibilities for generating gaugino mass nonuniversality. In ad- 
dition, there is a new feature in this string model, not present in GUT models, 
in that the universal scalar (mass) 2 at the unification scale, i.e. 7Bq, is depen- 
dent on 5gs which therefore correlates the universal scalar mass with the gaugino 
mass nonuniversality. An analysis of the maximum and the minimum curves for 
a xp (scalar) as a function of m x is given in Fig.6 when the parameters in the 
model are varied over a range of allowed values. A comparison with Fig.l shows 
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that a XP (scalar) can be larger than for the universal SUGRA case by a factor of 
as much as 10. One also finds that the range of the neutralino mass extends to 
about 575 GeV significantly beyond what one finds in the universal SUGRA case 
even without inclusion of the coannihilation effects. 

We discuss next dark matter for a class of D brane models. Over the recent past 
there has been considerable interest in the study of Type IIB orientifolds and their 
compactificationsjP], ^5], |3D|. We consider here models with compactifications on 



a six torus of the type T 6 = T 2 x T 2 x T 2 . In models of this type one has a 
set of 9 branes, 7* (i=l,2,3) branes, 5» branes and 3 branes. This set is further 
constrained by the requirement of N=l supersymmetry which requires that on has 
either 9 branes and 5j branes, or 7* branes and 3 branes. Model building on branes 
allows an additional flexibility in that one can associate different parts of the 



Standard Model gauge group with different branes. In Ref. [42, 46 a brane model 



using two five branes 5i and 52 was investigated while in Ref. [43 1 models using 



9 brane and 5j brane were investigated. We pursue here the implications of the 



latter possibility [|43|]. In one of the models of Ref.Q] the Standard Model gauge 
group is associated with the branes in the following way: the SU(3)c x U(l)y is 
associated with the 9 brane while SU(2)i is associated with the 5i brane. Further, 
it is assumed that the SU(2)r singlets are associated with the 9 brane while the 
SU(2)l doublets are associated with the intersection of 9 brane and 5i brane. The 
soft SUSY breaking sector of this 9 — 5i brane model is then given as follows: The 
gaugino masses rhi (i=l,2,3) corresponding to the gauge group 577(3), SU{2) and 
U(l) are parametrized by||43|| 

rhi = y/2>mz svn6e~ lls = rh 3 = —A Q , m 2 = V3wi3 cos#©ie~ 171 (37) 

while the SU(2)l singlets are parametrized by mg and SU (2)l doublets are parametrized 
by mg 5l where U4~3| 

m\ = m|(l -3 cos 2 eQj), m 2 g5i = m|(l - - cos 2 6(1 - 2 )). (38) 

Here the 9 and ©i are the directions of the Goldstino in the dilaton and the moduli 
VEV space as discussed earlier. To avoid generating tachyons in the theory one 
needs to impose the constraint cos 2 9Q\ < 1. A second version of this model was 
also discussed in Ref. [[43|[ . Here one associates SU(3) x U(l)y with the 5i brane 



while the SU{2) L with the 9 brane, and assumes that the SU{2) R singlets were 
associated with the 5i brane while the SU(2)l doublets are associated with the 
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intersection of 9 brane and 5i brane as before. The soft SUSY breaking sector 
of this model can be gotten from the model discussed above by the interchange 

cos#Oi < > sin 9. In the analysis of this paper we focus on the version of the 

model given by Eqs.37 and 38. As usual we assume that the parameter \i is free 
and we determine it via the constraint of the radiative breaking of the electro-weak 
symmetry. Again as in other cases we have considered we set the CP phases to 
zero. Eqs.37 and 38 show that one has nonuniversality at the unification scale both 
in the scalar sector as well as in the gaugino sector. The limit of universal scalar 
mass corresponds to 0i = l/\/3. Our focus in this paper is on the nonuniversality 
in the gaugino sector, and so for the numerical analysis we set Qi = l/v3. In this 
case the numerical analysis shows that the allowed neutralino mass range extends 
up to about 650 GeV. 

5 Gaugino mass nonuniversality and LSP mass 
range 

In SUGRA models with universal boundary conditions at the GUT scale, the 
allowed range of the LSP typically does not exceed 200 GeV and with imposition 
of additional constraints it is often significantly less. As discussed in Sees. 3 and 
4 in the presence of gaugino mass nonuniversality one finds that the allowed LSP 
range is increased. For the 1 + 24 case the allowed LSP range extends to about 
220 GeV for the case C24 = —0.1. For the 1 + 75 case one finds that C75 = —0.06 
gives an LSP range which extends to 250 GeV, while for the 1 + 200 case with 
C200 = 0.1 the allowed LSP range extends to 375 GeV. A similar situation holds for 
the string/brane models. Here one finds that the allowed LSP range can extend up 
to 600 GeV. These extended regions arise even without inclusion of coannihilation 



effects which is known to extend the allowed regions also up to about 700 GeVp6 
In the regions of the parameters space considered the effects of coannihilations 
would infact be negligible since we are considering only those configurations for 
which Aj > 0.1. The specific mechanism by which the neutralino mass range is 
extended is also different than in the case of coannihilation. Thus for the case of 
coannihilation the increase in the allowed LSP range occurs as a consequence of a 
coupled channel effect while in the the case of nonuniversalities the extension of 
the allowed region of the LSP arises due to a significant increase in the value of J 
for certain ranges of nonuniversalities. Thus one may expand J as the sum over the 
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final state channels in the XX annihilation so that J-J(ff) + J(WW) + J(ZZ) + 
J(Zh°)+J(ZH°) + etc. The region of the nonuniversality parameter space which 
leads to an enhancement of a xv can also lead to an enhancement of the cross- 
sections for the XX annihilation into the final states W + W~ , ZZ, Zh, ZH etc and 
hence to an increase in J which leads to a decrease in the relic density down 
to permissible limits consistent with constraints. Thus regions of the parameter 
space which would otherwise be excluded are now included when the gaugino mass 
nonuniversalities are included. 



6 Conclusion 

In this paper we have analyzed the effects of the nonuniversality of the gaugino 
masses on neutralino dark matter in SUGRA, string and D brane models under 
the constraint of R parity invariance. It is found that nonuniversality effects can 
enhance the x~P cross-section for scalar interactions by as much as a factor of 10. 
We also carried out an analysis of the uncertainties in the numerical determination 
of ff (i=u,d,s) and find that with the current state of uncertainties the x~P cross- 
section cannot be pinned down to better than a factor of about 5. Our analysis of 
the gaugino mass nonuniversality also exhibits another important phenomena, i.e., 
that the allowed range of the neutralino mass can be extended up to about 600 GeV 
even without inclusion of the coannihilation effects. The effects of coannihilation 
were not considered and inclusion of these effects may further increase the allowed 
neutralino mass range. The extended LSP mass range should be of interest in 



the experimental searches for dark matter in the current dark matter detectors [pO 
|31| , |47j and in the design of new dark matter detectors which are at the planning 
stagepH- 
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Appendix A: Effects of gaugino mass nonuniversality on sparticle masses 

In this appendix we give analytic solutions to the one loop renormalization group 
equations including the effect of gaugino mass nonuniversality. The analytic for- 
mulae for these with universal boundary conditions were given in Ref.[[49[ and for 



nonuniversality in the scalar mass sector in Ref. ||19||. Here we limit ourselves to 
the gaugino mass nonuniversality. These formulae are found useful in gaining an 
analytic understanding of the nonuniversality effects. The one loop RG formulae 
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are given in several papers (see, e.g. Ref . [|49|1 ) and we do not reproduce them here. 



Rather we discuss the solutions under the boundary conditions where is the 
universal scalar mass, rhi(0) (i=l,2,3) are the gaugino masses for the gauge group 
sectors U(l), SU(2), SU(3) and A is the universal trilinear coupling all taken at 
the unification scale. The Higgs mass parameters, the trilinear couplings, and the 
squark and slepton masses at the electro-weak scale are all sensitive to the effect 
of gaugino mass nonuniversality. The simplest case is that of the mass parameter 
m 2 Hi for the Hi Higgs which couples to the down quark. Here one finds 

m 2 Hl =m 2 + 5 (§/ a (t) + ^Ut))m\ (39) 

where 

^'^-tttW"— ><— > 2 < 40 » 

Pi (1 + (Jit) 1 a G mi 

Here t = ln(M^/Q 2 ), (3, = (6^/4^)^(0) where b { = (33/5, 1, -3) for 17(1), SU(2), 
and SU(3), and «i(0) = ai/An. The /j(i) contain the nonuniversality effects. The 
evolution of the mass parameter for the Higgs H 2 involves the evolution of the 
trilinear coupling in the stop channel at the electro-weak scale and this coupling 
is given by 

At{t ) = A ° + mi(H 2 - dY ° Hs ) (41) 
V ; 1 + 6F 7 ^ 1 + 6Y F J V ; 

where 

16- ~ 13- rt 
H 2 = ao( —h 3 + 3h 2 + —h,), H 3 = / E(t)H 2 (t) 

k = « (42) 
1 + (Jit a G mi 

2 

where Yq is the top Yukawa coupling at the GUT scale and the functions E and 
F are as defined in Ref.|49[]. We note that the first term on the right hand side of 
Eq.(40) which arises purely from the top Yukawa coupling evolution is unaffected 
by nonuniversality while the second term in affected through the modification of 
hi. For the mass parameter m# 2 for Higgs H 2 that couples with the top one finds 

m 2 H2 = m \e{t) + A m mif(t) + m 2 (h(t) - k(t)A 2 ) (43) 
where the functions h(t) and k(t) are unaffected by nonuniversality and are as 



given in Ref.|E| while the functions e and / are modified due to nonuniversality. 
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and 



e and / are given by 

6(t) - 2 [ D(t) + 3D(t)2 + Hs] (44) 

/(*) = -^T^; = (1 + 6W)) (45) 

Here the various tilde functions containing the nonuniversality are defined below 

G l {t)=F 2 {t)- l -{H 2 f 

G 2 (t) = 6F 3 (t) - hit) - AH 2 (t)H 4 (t) + 2F(H 2 ) 2 - 2H 6 (t) 

F 3 (t) = F(t)F 2 (t) - f dt'E(t')F 2 (t'), F 4 (t)= f dt'E(t')H 5 (t') 
Jo Jo 

H 4 (t) = F(t)H 2 (t) - H 3 (t), H 5 (t) = a G (-fh(t) + 6/ 2 (f) - 

H 6 (t) = f o dt'(H 2 ) 2 E(t'), H 8 (t) = &G(~Mt) + hH) - 1/0 (46) 

The squark and slepton masses are also affected by nonuniversality. For the up 
squarks in the first two generations one finds 



mljt) = m 2 + m 2 ui + a G [^f 3 + -f 2 + ^fi]m\ + {\~\ sin 2 9 W )M 2 Z cos 2(3 



8 8 2 
ml iR {t) =m 2 + m 2 ui + a G [-f 3 + j^fi]m\ + - sin 2 6 W M 2 cos 2(3 (47) 

The analysis of the first two generations of the down quarks and of the sleptons is 
similar. Finally we discuss the third generation squarks. Here one finds 



m~ L (t) = l -m\ + m 2 + ^mfj + +a G [p 3 + ^r/i]™| + (-\ + \ s ™ 2 W )M 2 cos 2(3 

8 2 1 

m~ (t) = m„ + Bij + 5;g[o/3 + Tr/i] m i — o sm2 @wM z cos 2/3 



1 2 

mf^ (f ) = raj + m t 2 + (- - - sin 2 6 W )M 2 Z cos 2(3 

2 

mf (t) = ml + m 2 + - sin 2 W M§ cos 2/5(48) 
fi 3 



where mjy and tjiq are defined by 



m u = \ m o + - jj^o + g/wng + [-e + o:g(-/3 - /2 + l/i)] 



™?q = \™l + \fAomi - \kA 2 + \hml + [\e + a G (^h +h~ ^r/i)] (49) 
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All the sparticle mass relations limit to the universal case j|9| when we set 
aij(0)mj(0) / aami = 1. 

Appendix B : nonuniversality effects on /x 

Since /x is determined via the constraint of electro-weak symmetry breaking it is 
sensitive to the gaugino mass nonuniversality. Thus one has 

^ = ( m 2 Hi - m 2 H2 tan (3 2 ) (tan (3 2 - l)" 1 - -Mf + A/x 2 (50) 
where A/x 2 is the loop correction. From the above one finds 

^- = (t 2 - l)-\m\g' r - t 2 (mle' r + A m Q m,f' r )) + ^ (51) 

UC r 2 2 2 OC r 

Here g' r = where g = a G (f/2 + and f r = J£. For large tan/3 Eq.(51) 

reduces down to 

dfi 2 , , <9A/x 2 

— — = -{m\e' r + A m mifr) + — (52) 

OC r 2 2 c/ C? , 

These results lead to Eq.(34). 



Tables: 



Table 1: Nonuniversalities at Mx- 



SU(5) rep 


n\ 


n r 2 




1 


1 


1 


1 


24 


-1 


-3 


2 


75 


-5 


3 


1 


200 


10 


2 


1 



Table 2: Uncertainties in a n N, do and x 



Ref. 


ow )JV (MeV) 


Ref. 


ct (MeV) 


Ref. 


X 




a 


45 ±8 





35 ±5 


O'o/^ttN 


0.74 ±0.25 






48 ± 10 


H 


36 ±7 


H 




0.64 ±0.03 




37| 




50 ± 10 


^0 


35.5 ±6 


M 




0.57 ±0.1 




28j 




49 ±8 






X 


0.67 ±0.18 




48 ±9 











Figure captions 

Fig.l: Exhibition of the dependence of a xp (scalar) on the uncertainties in f%, f v d 
and fl for the universal case. The plots are for three sets of ff discussed in Sec. 3. 
The parameter space spanned is as discussed in the text. 
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Fig. 2: Plot of gaugino-Higgsino components X h q (n=l,2,3,4) as a function of c 2 4 
for the input data m = 51 GeV, m x =70 GeV, tan/3 = 10, A t /m = —7 where 
X20 = —|X2o I and X40 = — 1-^40 1 • 

Fig. 3: Plot of the maximum and the minimum curves a xp (scalar) as a function 
of m x for the case when one considers admixtures of 1 + 24 representations with 
c 2 4 = —0.1 (dashed), c 2 4 = (solid) and c 2 4 = 0.08 (dotted) when the other pa- 
rameters are varied over the assumed naturalness range as discussed in the text. 

Fig. 4: Plot of the maximum and the minimum curves a xp {scalar) as a function 
of m x for the case when one considers admixtures of 1 + 75 representations with 
c 75 = —0.06 (dashed), C75 = (solid) and C75 = 0.04 (dotted) when the other 
parameters are varied over the assumed naturalness range as discussed in the text. 

Fig. 5: Plot of the maximum and the minimum curves a xp (scalar) as a function of 
m x for the case when one considers admixtures of 1 + 200 representations with 
C200 = 0.1 (dashed), c 2 oo = (solid) and c 2 oo = —0.08 (dotted) when the other 
parameters are varied over the assumed naturalness range as discussed in the text. 

Fig. 6: Plot of the maximum and the minimum curves for a xp (scalar) as a function 
of m x for the heterotic string model O-II. The range of parameters consists of e, e' 
in the range 0.0025-0.01, 9 in the range 0.1-1.6, Sqs hi the range -1 to -10, m 3 / 2 in 
the range up to 2 TeV, and values of tan/3 range up to 25. 

Fig. 7: Plot of the maximum and the minimum curves for a xp (scalar) as a function 
of m x for the 9 — 5i D brane model. The range of parameters consists of 9 in the 
range 0.1 — 1.6, m 3 / 2 in the range up to 2 TeV, and values of tan j3 range up to 25. 
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